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)0 ABSTRACT 

o. 

\0 . We analyze the structure and stability of the transition layer (or front) that connects 

the cold neutral medium and warm neutral medium in the plane-parallel geometry. 
p ! Such fronts appear in recent numerical simulations of a thermally bistable interstel- 

lar medium. The front becomes an evaporation or condensation front depending on 
the surrounding pressure. The stability analysis is performed in both long- and short- 



O 

H 



c^ ! wavelength approximations. We find that the plane-parallel evaporation front is un- 

stable under corrugational deformations, whereas the condensation front seems to be 
k> , Stable. The instability is analogous to the Darrieus-Landau instability in combustion 

H I front. The growth rate of the instability is proportional to the speed of the evapora- 

tion flow and the corrugation wavenumber for modes with wavelength much longer 
than the thickness of the front, and it is suppressed at scales approximately equal to 
the thickness of the front. The timescale of the instability is smaller than the cooling 
timescale of the warm neutral medium (~ 1 Myr), and can be as small as the cooling 
timescale of the cold neutral medium (~ 0.01 — 0.1 Myr). Thus, this instability should 
be one of the processes for driving the interstellar turbulence. 

Subject headings: headings: hydrodynamics — ISM: kinematics and dynamics — 
method: analytical — instabilities 



1. INTRODUCTION 

It is widely known that the low- and mid-temperature parts of the interstellar medium (ISM) 
are composed of a thermally bistable fluid that results from the balance of radiative cooling and 
heating due to external radiation fields and cosmic rays (Field, Goldsmith & Habing 1969; Wolfire 
et al. 1995, 2003). The two stable phases are called the cold neutral medium (CNM) with tem- 
perature T^IO — lO^K and the warm neutral medium (WNM) with temperature T ^ 10^ K. 
Zel'dovich & Pikel'ner (1969) and Penston & Brown (1970) studied the transition layer or front 
that connects these stable phases assuming plane-parallel geometry. They found that a static so- 
lution is obtained at the so-called the saturation pressure, and steady phase exchange (evaporation 
and condensation) solutions are obtained when the pressure deviates from this value. 

From an observational point of view, interstellar clouds are characterized by suprathermal 
emission-line widths, which can be related to the supersonic internal motion of gas, or "turbu- 
lence," in clouds. Recent numerical simulations have shown that supersonic turbulence in an 
isothermal or an adiabatic medium decays quickly, that is, within a crossing time, irrespective 
of the effect of magnetic fields (Mac Low et al. 1998; Stone et al. 1998; Ostriker et al. 1999). In 
contrast, Koyama & Inutsuka (2002) found long-lasting turbulent motion of tiny cloudlets in the 
WNM in a numerical simulation of the two-phase medium, which are formed as a consequence 
of thermal instability (see also Koyama & Inutsuka 2004, Inutsuka & Koyama 2004). Audit & 
Hennebelle (2005) found similar energetic turbulent motions in colliding flows, a topic has been 
studied further by Heitsch et al. (2005) and Vazquez-Semadeni et al. (2006). However, a detailed 
analysis of the mechanism that maintains such turbulent motions remains to be done. 

The two-phase evaporation front is analogous to a combustion front. The latter are unstable 
to corrugational deformations, a phenomenon known as the Darrieus-Landau instability (DLI; see, 
e.g.. Landau & Lifshitz 1987; Zel'dovich et al. 1985). From experiments and numerical simula- 
tions of chemical and nuclear flames, the DLI is known to develop into turbulence in the nonlinear 
regime, and this turbulence leads to accelerated combustion (Sivashinsky 1983; Niemeyer & Hille- 
brandt 1995; Blinnikov et al. 1995). It is also known that there is a similar instability that applies 
to evaporation fronts under the discontinuous-front approximation (Aranson et al. 1995; Inutsuka 
et al. 2005). Thus, the instability is expected to be a driving mechanism of turbulence in the 
two-phase medium. 

However, in the framework of the discontinuous front approximation, we cannot determine the 
most unstable scale of the instability and its growth timescale owing to the neglect of the thickness 
of the front. Therefore, in this paper we undertake a linear stability analysis of the two-phase front 
by accounting for the effect of the finite thickness of the front. 

The paper is organized as follows. The steady smooth fronts that describe the structure of the 



transition layers are calculated numerically in §2. In §3, we briefly review the linear stability of the 
front under corrugational deformations in the framework of the discontinuous-front approximation, 
which corresponds to analysis in the long-wavelength limit. Then the linear stability of the smooth 
transition layers that are obtained in §2 is studied under the isobaric approximation which gives 
an appropriate result for small wavelength modes. Finally, we summarize our results and discuss 
their implications in §4. 



2. STRUCTURE OF STEADY FRONTS 

We start the with fluid equations describing the dynamics of an unmagnetized, optically thin 
ideal gas of density p, temperature T, pressure p, and velocity v that is heated externally and cooled 
radiatively, 

^rr + ■^^^;—{p^i^k + p^ik) = , (2) 

Ot OXk 

^^ ^ \e + p)v,-k^\ = -pC{p,T), (3) 



dt dxk { dxk , 
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P=-pT, (4) 

i^-^ + ^, (5) 
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Here C is the net cooling function, k is the thermal conductivity, 7 is the ratio of specific-heats, R 
is the gas constant, and fi is the mean molecular weight. We adopt the following simplified cooling 
function 

pC = n{-r + nA), (6) 

where 

r = 2 X 10~26 ergs s~^ , (7) 

,1 / -118400 \ „ /-92\ o 1 

A = 7.3 X 10-21 exp(^^^-^^j+ 7.9 X 10-2^ exp(^^j ergs cm^s-^ (8) 

Figure 1 shows the equilibrium state of this cooling function (£(n, T) = 0) in the number density- 
pressure plane. A two-phase structure is possible in the case that lO^'^ < p/ks ^ lO'^'i K cm"'^ 
for this simplified choice of the cooling function. 



In this section, we consider the steady state transition layer (or front) that connects the CNM 
and WNM in plane-parallel geometry, where the CNM and WNM are in a state of radiative equi- 
librium {C = 0). From equations (l)-(3), the state on either side of the front must satisfy the 
following conservation laws: 

[pv,] = 0, (9) 

[p + pvl]=0, (10) 

[iE + p)v^]=Q, (11) 

where brackets represent a difference across the front {[F] = F(x = oo) — F{x = — oo)) and 

Q = - f pCdx, (12) 

is the net energy acquired as a result of external heating and radiative cooling. Shchekinov & 
Ibanez (2001) give a detailed study of the solutions of equations (9)-(l 1). In the ISM, the following 
three types of solutions are important: 

1. Evaporation.-The CNM acquires energy through the front and evaporates to the WNM. This is 

obtained by choosing the left- and right-hand states to be the CNM and WNM, respectively, 
with Q>0. 

2. Condensation.-The WNM loses energy through the front and condenses to the CNM. This is 

obtained by choosing the left- and right-hand states to be the WNM and CNM, respectively, 
with Q <0. 

3. Saturation.-This is the static solution for which external heating and radiative cooling are bal- 

anced inside the front {Q = 0). 



2.1. Solutions Connecting Two Radiative Equilibrium States 

In this section, we calculate the structure of the fronts that connect the CNM and WNM. The 
energy equation (3) can be rewritten into the form 

-i R dT dp , . 



7 — 1 yU dt dt 

where d/dt = d/dt + v ■ V. We seek a one-dimensional solution that satisfies thermal equilibrium 
at infinity: 

£(pi,Ti)=0, (14) 

>C(P2,T2) = 0, (15) 



where subscripts 1 and 2 denote the values at x = — cxd and x = oo, respectively. Omitting the 
y— , z— and time dependences of the variables, and integrating equations (1) and (2) with respect 
to X, we obtain the basic equations corresponding to the three conservation laws: 
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where v is the x-component of the velocity. From equations (16), (17), and (4), the velocity, pres- 
sure, and density can be expressed in terms of the temperature, mass flux j and total momentum 

M: 
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The density pi and temperature Ti at x = — oo are obtained from the solution of equation (14) if 
the pressure at x = — oo (pi) is given. Then the total momentum can be expressed as a function of 
the mass flux and given pressure at x = — oo, M = M{j, pi), and the boundary conditions 
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(22) 
(23) 



are imposed on equation (18). Equation (18) is a second-order ordinary differential equation with 
respect to T, whereas three boundary conditions (22) and (23) are imposed on it. Therefore, it has 
a solution only when the mass flux j is tuned as an eigenvalue of the solution, and j is expressed 
as a function of pi (Shchekinov & Ibanez 2001). 



In the case of a static front, that is, j 

up C dT 



0, equation (18) can be transformed as follows: 
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(Zel'dovich & Pikel'ner 1969, Penston & Brown 1970). Equation (24) expresses the balance be- 
tween cooling and heating inside a static front. The pressure that allows such a static structure 



is called the saturation pressure. Substituting equations (6)-(8) into equation (24) and solving, 
we obtain a saturation pressure Ps/ks = 2612 K cm~^ for our choice of cooling function. If pi 
is higher (lower) than this pressure, fluid elements that pass through the front experience a net 
cooling (heating) inside the front. Thus, such a front becomes a condensation (evaporation) front. 

We numerically solve equations (16)-(18) with the boundary conditions of equations (22) and 
(23). We take the CNM to be state 1. The dominant contribution to thermal conductivity is due 
to neutral atoms, so we use k, = 2.5 x 10'^ T^/^ ergs s"^K~^cm~^ (Parker 1953). We plot the 
solutions in the number density-pressure diagram and in terms of temperature structure in Figures 
1 and 2, respectively, for values pi/k-Q of 1500, 2000, 2612, and 4000 K cm^'^. The structure of 
the temperature, pressure, number density, and velocity for pi/k-Q = 1500 K cm"^ are plotted in 
Figure 3. 

Figures 1 and 3 show that the structure is almost isobaric. That is due to the small Mach 
number of the flow (typically M. — 10"^). From equation (16) and (17), the pressure difference 
between the CNM and WNM is 

Pav Pw 

~ 7 (— ) =10-^ (26) 

where subscripts w and c denote the values in the WNM and CNM, respectively, c = (7 i? T/yu)^/^ 
is the sound speed, and we use the fact that v^ ^ Vc. 

From equation (18), the thickness of the front is determined by the Field length Ip, 

kT 
Axr^ J—^=If (27) 



(Begelman & McKee 1990; Ferrara & Shchekivov 1993), which is a function of local variables 
and therefore takes different values in the CNM and WNM: 

If,c = V^TjnlA ~ 10-2 pc , (28) 

If,^ = ^KT^/n^T ~ 10-^ pc , (29) 

This property can be seen in Figures 2 and 3. The total thickness is essentially determined by the 
Field length in the WNM, If,w 

In Figure 4, we plot the mass flux j and the velocity of the WNM t;„ as functions of pi in the 
top and bottom panels, respectively. The typical evaporation flow speed is on the order of 0.1 km 
s"^ for this type of solution connecting two radiative equilibrium states. 



2.2. Solutions with Finite Spatial Extent 

We can obtain a larger evaporation flow speed (or smaUer condensation flow speed), if we 
change the boundary condition in equation (22) to a higher CNM temperature, without changing 
Pi, and impose the CMN boundary conditions at a finite distance from the front. When such 
conditions are imposed, the CNM is no longer in a state of thermal equilibrium: £(pc, ^c) 7^ 0. 
The top and bottom left panels of Figure 5 show, respectively, the temperature structure and n — p 
diagram for solutions of infinite and finite spatial extent with Pc/ks = 1200 K cm~'^. The CNM 
temperature and evaporation flow speed are Tc = 250K, and v^ = 0.34 km s^^ for the finite spatial 
extent solution and Tc = 35. 5K and v^ = 0.14 km s~^ for infinite one. In both solutions, the WNM 
temperatures are the same. From the bottom left panel of Figure 5, one can see that the net heating 
of a fluid element that passes through the front is stronger for the finite spatial extent solution 
than for the infinite one, because the cooling becomes weaker. Thus the evaporation flow speed 
becomes faster when the CNM temperature is larger than the thermal equilibrium temperature. 

Of course, we can obtain a lower WNM temperature solution if we impose the boundary 
conditions for the WNM at a finite distance from the front. In this case, the condensation flow 
speed becomes large (or the evaporation flow speed becomes small), because of the stronger net 
cooling (weaker heating) inside the front. The top and bottom right panels of Figure 5 show, 
respectively, the temperature structure and n — p diagram of the infinite and finite spatial extent 
solutions for Pc/k-Q = 4000 K cm~'^. The WNM temperature and condensation flow speed are 
r„ = 5000K, and v„ = — 3.73 x 10"^ km s~^ for the finite spatial extent solution and T^ = 81 15K, 
and fw = —3.02 x 10"^ km s~^ for the infinite one. In both solutions, the CNM temperatures are 
the same. These kinds of solutions are naturally realized in dynamical environments. 

We plot the evaporation flow speed t>„ versus the temperature of the CNM in the case that 
Pi/ks = 1000, 1200, 1500, and 2000 K cm~^ in Figure 6. The evaporation flow can reach speeds 
as large as 1 km s~^ In §3, we show that the growth rate of the DLI (eq. (48)) is proportional to 
the evaporation flow speed v^. Thus, this effect makes the DLI more effective. 



3. STABILITY OF PHASE TRANSITION LAYERS 

3.1. Long- Wavelength Approximation 

In this section, we focus on perturbations on the scales much larger than the thickness of the 
front (the Field length), so we treat the transition layer as a discontinuous front. The effect of the 
finite thickness of the transition layer is considered in the §3.2. 

We assume that the discontinuous front is located in the plane x = and that the direction 



of the flow is parallel to the x-axis, and we choose a reference frame that moves together with 
the front. The unperturbed state satisfies the jump conditions (9)-(ll). As stated in the previous 
section, the jump conditions admit three types of solution, that is, evaporation, condensation, and 
saturation solutions. In the following, we carry out a linear stability analysis of these fronts under 
a corrugational deformation. The situation is illustrated schematically in Figure 7. We assume that 
perturbations of the fluid on both sides of the front are incompressible, because the flow speed is 
sub-sonic in either region. The basic perturbed equations of the fluid are 

kx Svx + ky 6vy = , (30) 

pi,2 (^ Svx + kx vi^2 Svx) = -k^5p, (31) 

Pi^2{^Svy + kr,vi^2Svy) = -ky6p, (32) 

where subscript 1 and 2 denote the left (x < 0) and right (x > 0) sides of the unperturbed state, 
respectively, and we assume that the perturbation is proportional to e* ^'^'^ ^^^^ ^'"^ *^ . From equations 
(30)-(32), the perturbations in the pressure and the a;-component of the velocity can be expressed 
in terms of the y-component velocity perturbation: 

. pi {ki^^ vi-u) (_) 

5p = 6v^y > , (33) 

Ky 

5vx = -^Sv^-\ (34) 

Kx 

6vy = 5^-) , (35) 



on the left side of the front, and 

k 



. P2 {kx V2-Uj) (,^ P2 {kx V2 - Uj) (^. 

6p = -^-^ 5vl;^> - ^-^ — '- 6vlf> , (36) 

y y 

^^. = -^H*'-i*^i", (37) 

thx rhx 

6vy = 54+) + 6vl;^ , (38) 

on the right side, where kx = —i ky, kx = i ky, and kx = w/vr. Each of these modes damps 
as it propagates downstream (the (+) and the (s) modes) or upstream (the (— ) mode), where we 
have classified the (s)-mode as a downstream component because we seek an unstable solution 

(lm[u] > 0). 

Let us consider the matching conditions for the perturbations at the front. The perturbation in 
the position of the front can be expressed as 5xf oc e* ^''^ y-^^^) ^ and thus the velocity perturbations 
normal and tangential to the front are given by 

5vn = 5vx + ioj 5xf , (39) 

5vt = Svy + ikyV 6xf . (40) 



The matching conditions of the front are perturbed version of mass and momentum flux conserva- 
tion laws: 



[p5Vn] = 0, 

[6p + 2pv6vn] =0, 
[p V 6vt] = . 



(41) 
(42) 
(43) 



We modify the equation (41) for the following physical reasons. The discontinuous-front approx- 
imation can be justified when we consider a perturbation whose wavelength is longer than the 
thickness of the front; the structure of the front is hardly deformed under such a perturbation. 
Thus, the mass flux passing through the front will not changed, because it is determined by the 
structure of the front (eq. (12)) 

Pi 5Vn,l = p2 SVn,2 = . (44) 

Substituting equations (33)-(40) into equations (42)-(44), we obtain the characteristic equation 



where 



V 



Mr7 = 0, 
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1 -1 -1 



\ 

i kyV2/uj 



(47) 



ijj 

1 

\ iky{rd, - 1)^2 

and r^ = P2I p\ = f i/f2 is the compression ratio between regions 1 and 2. Taking the determinant 
of M, we obtain the dispersion relation for the front. 



/ 



-ra + rdJl-Td + r^^ 

-lUJ = ^ V2 ky. 

l + rd 



(48) 



In the case of evaporation, r^ = Pw/Pc < 1, and in case of the condensation, r^ = pc/pw > 1- 
Thus, the evaporation front is unstable and the condensation front is stable. Note that this is a 
similar to the Darrieus-Landau instability in combustion fronts (see, e.g.. Landau & Lifshitz 1987). 
The similarity of the evaporation and the combustion fronts can be understood if one imagines the 
CNM as cold fuel, the WNM as hot exhaust, and Q as the energy released by combustion. 

From equation (48), the growth rate of the instability is proportional to the wavenumber of 
the corrugation ky. Therefore, the most unstable scale is infinitesimal, and we cannot estimate the 
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physical scale of the instability or its growth time. This due to our neglect of the thickness of the 
front. 

To obtain a more realistic dispersion relation, Aranson et al. (1995) considered the effect of 
the curvature of the front (see, Graham & Langer 1973; Aranson et al. 1995; Nagashima et al. 
2005) and changed the right-hand side of equation (44) to a ky, where a is a mass flux variance 
due to the curvature effect. This stabilizes the front for the short wavelength modes, and is also 
known as the "Markstein effect" in the case of a combustion front (Markstein 1964; Zel'dovich et 
al. 1985). 

However, the most reliable way to determine realistic dispersion relation is to consider the 
thickness (characteristic scale) of the front in the analysis. By doing so, we can account for the ef- 
fects of thermal conduction and thermal stability inside the front that are lost in the discontinuous- 
front approximation. In the following section, we therefore study the effect of the thickness of the 
front. 



3.2. Short- Wavelength Approximation 

In this section, we study the linear stability of the transition layer obtained in §2.1. However, 
it is difficult to carry out the analysis without approximations. Thus, we consider a perturbation 
whose wavelength is much shorter than the acoustic scale la, which is defined as c^ tc, where Cg is 
the sound speed and tc is the characteristic cooling (or heating) timescale. Normally, the acoustic 
scale is longer than the Field length: /„ ^ ~ 10 pc in the WNM and /„ c ~ 0.1 pc in the CNM. For 
such small-scale modes, pressure balance sets in rapidly, and perturbations of the thermodynamic 
variables (temperature and density) are important. On the other hand, perturbations of the fluid 
dynamical variables (pressure and velocity) are important, for the large-scale modes considered in 
§3.1. 

We assume an isobaric perturbation of the front. The basic equations are perturbed version of 
equations (13) and (4) for the themodynamic variables: 

^ ^{po{n6T + voD5T)+vo{DTo)Spo} = 



7-1 ^ 
6T 6p 



{DTo)^ , {D'To) {DT,^' 



k{To) ( {D' - kl) + ^^D + ^^^ - ^^^ ]dT-Lr6T- L, 6p , (49) 



^ , 0, (50) 

where a subscript zero denotes an unperturbed variable, D = d/dx is a differential operator, and 
we write Lt = d {p C)/dT and Lp = d (p C)/d p. We assume that a perturbed physical variable 



11 



/ takes the form 

f{r,t) = fo{x) + Sf{x) exp{ikyy + nt). (51) 

Since the unperturbed state is spatially uniform at infinity, we can obtain asymptotic solutions 
by assuming solutions proportional to e^^'. From equations (49) and (50), x is determined by the 
following polynomial: 

i^{Tc,w) ix^ - kl) - pc,w 7 — {n + w^c,w X) - ^t|c,w + ;7T^^p|c,w = . (52) 

We seek a localized solution around the front that vanishes at infinity. Therefore, the asymptotic 
boundary conditions are 

5T{x) = Cie^^% (53) 

D6T{x) = CiXce^^% (54) 

far from the front in the CNM, and 

6T{x) = C2e''-% (55) 

D6T{x) = C2Xwe^"^ (56) 

far from the front in the WNM, where Ci and C2 are constant, and Re[xc] > and Re[xw] < 0. 
To avoid the freedom to translate the amplitude of the perturbation, we fix Ci. A solution can be 
obtained if ky is given and n and C2 are eigenvalues of equations (49) and (50) with the above 
boundary conditions. 

This eigenvalue problem has an infinite number of solutions (eigenmodes). In the analysis, 
we are not interested in overstable modes. Therefore, we search for a solution whose growth rate 
n is real. We numerically solve the eigenvalue problem for a given wavenumber. The dispersion 
relations of the most unstable eigenmodes are plotted in Figure 8 for Pc/ks = 1500, 2000 , 3000 
K cm"'^. In the case of evaporation, there is a range of unstable wavenumbers, whereas all values 
are stable in the case of condensation. In all cases, the eigenfunctions do not have a node; that is, 
the solutions are fundamental modes. The timescale of the instability is roughly on the order of 
the cooling timescale of the WNM tc,w ~ 1 Myr or shorter. The instability is stabilized at a scale 
approximately equal to the thickness of the front, If,w ~ 0.1 pc, as a result of thermal conduction. 
In the isobaric approximation, results for wavelengths comparable to or longer than the acoustic 
scale are not valid. We expect that the growth rate would be suppressed at longer wavelengths as 
in equation (48), if we were to not use the isobaric approximation. 
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4. SUMMARY AND DISCUSSION 



We have performed a linear stability analysis of the steady transition layer that connects the 
CNM and WNM in both the long- and short-wavelength limits. Our results show that the evapora- 
tion front that forms when the surrounding pressure is lower than the saturation pressure is unstable 
under corrugational deformations, whereas condensation fronts seem to be stable. 

The result of the long-wavelength (discontinuous front) approximation, equation (48), shows 
that the growth timescale of the instability tgrow is inversely proportional to the speed of the evap- 
oration flow v„ and the corrugation wavenumber ky-. 

where A = 2 n/ky is the wavelength of the corrugation. A typical density contrast in the ISM is 
about a hundred, r^ = 0.01. The results of the short-wavelength (isobaric) approximation show 
that the instability stabilizes at the scale of the thickness of the front, ~ 0.1 pc, which is essentially 
determined by the Field length of the WNM. 

We plot the dispersion relation for both the long- and short-wavelength approximations at 
a value of Pc/ks = 1000 K cm~^ in Figure 9. The dotted and dashed lines are, respectively, 
the dispersion relations in the long- and short-wavelength approximations. The solid line is a 
schematic dispersion relation expected to be obtained from an analysis without approximations. 
From this figure, we expect the most unstable scale to be about twice the stabilized scale. For such 
a scale, the growth timescale of the instability is comparable to the cooling timescale of the WNM 
(~ 1 Myr). The growth timescale can be shorter in a dynamical environment because of the effect 
of finite spatial extent as studied in §2.2. The evaporation flow speed becomes as large as 1 km 
s~^. For this flow speed, from equation (57), the growth timescale can be as small as the cooling 
time-scale of the CNM (~ 0.1 — 0.01 Myr). Therefore, the instability is expected to grow within a 
dynamical timescale of the ISM (~ 1 — 10 Myr) and may drive a turbulent two-phase medium in 
the nonlinear stage, as in the case of a combustion front. 

Koyama & Inutsuka (2002) and Audit & Hennebelle (2005) have reported on the two-phase 
turbulence in numerical simulations of the thermally bistable medium, for which the velocity dis- 
persion is on the order of the observed velocity dispersions in the diffuse ISM. They propose that 
turbulence in the ISM can be driven in a thermally unstable gas constantly supplied by shock 
propagation (Koyama & Inutsuka 2002) or converging flows (Audit & Hennebelle 2005; see also 
Heitsch et al. 2005; Vazquez-Semadeni et al. 2006). These shocks or converging flows are ubiq- 
uitously generated, for example, by supernovae, and correspond to energy input on large scales. 
In contrast, our analysis of the stability of the transition layer shows another aspect of the dynam- 
ics of the multi phase medium that may be created by the thermal instability: the instability of 
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the evaporating front can generate fluctuating motions even without external mechanical forcing. 
Note, however, that the maximum growth rate of the instability depends sensitively on the velocity 
of the evaporating flow. A further detailed study of the nonlinear behavior of the instability will be 
presented in our next paper. 
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Fig. 1. — Equilibrium state with the cooling function given by eqs (6)-(8) (thick solid line). The 
thin solid line shows the steady solution of a saturation front (p/Zcb = 2612 K cm"'^), the evap- 
oration solutions for pi/k-Q = 1500, and 2000 K cm~^ are plotted as dotted and dashed lines, 
respectively. The condensation solution for pi/fce = 4000 K cm~^ is plotted as a dot-dashed line. 
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Fig. 2. — Temperature structure of the steady solutions. The solid, dotted, dashed and dot-dashed 
lines show the solutions for pi/k-Q = 1500 (evaporation), 2000 (evaporation), 2612 (saturation), 
and 4000 (condensation) K cm"'^, respectively. 
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Fig. 3. — Temperature T (solid line), pressure -p (dotted line), number density n (dashed line) 
and velocity v (dot-dashed line) structures of the steady evaporation solution for pi/Zce = 1500 K 
cm~^, where the number density has been multiplied by 100. 
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Fig. 4. — Mass flux j (top) and velocity of the WNM v^ (bottom) vs. pressure of the CNM pi. 
Positive velocities correspond to evaporation. 
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Fig. 5. — Temperature structure (top) and n — p diagram (bottom) for the finite spatial extent 
solutions. The left panels show the evaporation solution for p/ks = 1200 K cm~'^, and the right 
panels show the condensation solution at which p/ks = 4000 K cm^'^. We also plot the infinite 
spatial extent solutions connecting the two radiative equilibrium states (dotted lines). 



6 

W 0.5 



■TJ 0.4 

a; 
a; 
ft 

^ 0.3 

o 

.2 0.2 

'o 

ft 

03 0.1 

> 

a; 



Pj/k]3 - 1000 □ 

1200 o 

1500 A 

2000 o 



>oooo«o' 



«oooO 



ooO» 



oO» 



o.o-'^ 



temperature of CNM : Tc [ K j 



Fig. 6. — Evaporation flow speed v„ vs. temperature of the CNM in the case that Pc/k-Q = 1000 
(squares), 1200 (circles), 1500 (triangles), and 2000 (diamonds) K cm""'. Filled points correspond 
to the thermal equilibrium temperature of the CNM. 
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Fig. 7. — Schematic illustration of the perturbed front. We denote the wavenumber of the per- 
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Fig. 8. — Dispersion relations for the isobaric perturbation. The solid, dotted, dashed, dot-dashed, 
and dot-dot-dashed lines correspond to the unperturbed states of Pc/k^ = 1000 (evaporation), 
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respectively. 
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Fig. 9. — Dispersion relations of both the long- (dotted line) and short- (dashed line) wavelength 
approximations for j9c/A;b = 1000 K cm~^. The unperturbed flow speed is f„ = 0.23 km s~^, and 
the density contrast of the CNM and WNM is r^ = 0.042. The solid line is a schematic dispersion 
relation that is expected to be obtained from the analysis without approximations. 



